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Abstract 

We investigate the chaotic inflationary model using the two-loop effective potential of a self¬ 
interacting scalar field theory in curved spacetime. We use the potential which contains a non-minimal 
scalar curvature coupling and a quartic scalar self-interaction and which was found in Ref. [l|. We 
analyze the Lyapunov stability of de Sitter solution and show the stability bound. Calculating the 
inflationary parameters, we systematically explore the spectral index n s and the tensor-to-scalar ratio 
r, with varying the four parameters, the scalar-curvature coupling £o, the scalar quartic coupling Ao, 
the renormalization scale /r and the e-folding number N. It is found that the two-loop correction on 
n a is much larger than the leading-log correction, which has previously been studied in Ref. [2\ We 
show that the model is consistent with the observation by Planck with WMAP Em and a recent 
joint analysis of BICEP2 [5], 
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1 Introduction 

The measurements of the cosmic microwave background (CMB) fluctuations become increasingly impor¬ 
tant from the perspective of not only cosmology but also elementary particle physics. Useful indicators 
of the CMB fluctuations are given by the scalar (or density) fluctuations, S, the spectral index, n s , and 
the tensor-to-scalar ratio, r. In the inflation scenario a non-vanishing potential energy density of a scalar 
field induces exponential expansion on the universe. The origin of the CMB fluctuations is found in the 
quantum fluctuations of the scalar field. It is expected that these inflationary parameters restrict the 
models of particle physics. Although we have observed only one elementary scalar field, i.e. Higgs, it is 
quite natural to assume that other scalar fields exist and play a decisive role for the energy density and 
its fluctuation at early universe. 

In this paper we consider that the inflaton field is a real scalar field with a quartic self-interaction and 
a non-minimal scalar-curvature interaction at high-energy scale and study a possible model consistent 
with the CMB fluctuations at the two loop level. Inflation is thought to occur near the Planck scale. 
In such high energy scale the quantum correction may have some remarkable effect on the inflationary 
parameters. It is known that the spectral index and the tensor-to-scalar ratio are independent of the 
scalar quartic coupling, Ao, at the tree level. It is also known that there is an attractor on the (n s ,r) 
plane. The inflationary parameters, n s and r, converge to their universal model-independent values at 
the large scalar-curvature coupling limit El 13- Note that RG behavior of scalar curvature coupling £ is 
defined by the behavior of the corresponding quantum field theory at high energy (see: mm) so that it 
maybe tend to large or small asymptotic value at high-energy limit. 
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The inflationary parameters have been investigated up to the leading log level with respect to the 
scalar quartic coupling in Ref. [2]. The quantum corrections introduce the quartic coupling dependence 
for n s and r, but do not alter the attractor behavior. The standard model (SM) Higgs inflation has been 
investigated up to the next to leading log level in Refs. mm- In SM the scalar quartic coupling is 
extremely suppressed near the Planck scale. It has been pointed out that a large non-minimal scalar 
curvature coupling is necessary to reproduce the observed Higgs mass, n s and r. The remark is in order. 
Quantum field theory in curved spacetime induces log terms in the scalar four-point as well as in the 
non-minimal scalar-curvature sector, for the corresponding effective potential mug. The account of 
such quantum-corrected terms in the potential, especially RG improved effective potential, is done for the 
study of inflation in [TU HU HU HU HU HU HU HU HU HU HU HU HU • It is also interesting to note that the 
reconstruction of the inflationary scalar potential as is done in m maybe applied to such non-minimal 
inflationary scalar potential. 

The paper is organized as follows. Following Ref.jl. we introduce the effective Lagrangian up to the 
two-loop level in Sec. 2. The Lagrangian in the Einstein frame is formulated by the Weyl transformation. 
We show the concrete expressions for the scalar fluctuations, the spectral index and the tensor-to-scalar 
ratio at the two loop level in Sec.3. The Lyapunov stability of de Sitter solutions is studied in Sec. 4 in 
close analogy with the corresponding study for log-corrected higher-derivative quantum gravity [28]. In 
Sec.5, we numerically evaluate the inflationary parameters and show possible parameters consistent with 
the observed data. Finally we give some concluding remarks. 


2 Two-loop effective Lagrangian in Einstein frame 


Here we consider a massless scalar field with a non-minimal scalar-curvature coupling, £oj and adopt a 
simple chaotic inflation scenario near the Planck scale. We start from a Lagrangian density, 


£ (J) = + , 
y( J ) = 


% 4 , 

24^ ’ 


( 2 . 1 ) 

( 2 . 2 ) 


where the superscript (J) denotes the Jordan frame, g is the determinant of the metric tensor, g M „, £o> 
and Ao represents the scalar-curvature and scalar quartic couplings, respectively. The Jordan frame is 
characterized by the existence of the £-term. Here the reduced Planck mass is set as M p = ( 87 rG ) -1 / 2 = 1. 

In Ref. p] the closed expression for the two-loop effective potential is given under the linear curvature 
approximation, 
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where g represents the renormalization scale. It should be noted that some terms in the effective potential 
vanish if we set the scalar-curvature coupling as £0 = 1/6, called the conformal one. From Eas. (l2.1|) - (12.3I1 
we define effective scalar-dependent couplings, £ and A, as 
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The Lagrangian, Cpp, has a similar form with (12.11) thanks to the definitions (12.41) and (12.51) . This 
potential has the renormalization scale dependence which stems from the radiative corrections for <^ 4 
theory in curved spacetime. 

For calculations of the inflationary parameters, it is more convenient to change the frame into the 
Einstein frame where the £-term disappears. In order to change the frame, we consider the Weyl trans¬ 
formation, 

r v = fr 2 w, ( 2 . 8 ) 

where g is the metric tensor in the transformed frame. The Weyl factor, 12, is an analytic function with 
respect to the space-time coordinate. After this Weyl transformation the two-loop effective Lagrangian 
is rewritten as 
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where g 7 R and □ are the determinant of the metric tensor, the Ricci scalar and the d’Alembert operator 
in the transformed frame, respectively. We can transform the Jordan frame into the Einstein frame by 
choosing the Weyl factor in order that the £-term is eliminated. The suitable choice is 

12 2 = 1 + ^ 2 . (2.10) 

Then we redefine the scalar held to obtain the canonical kinetic term for the scalar held. The redefined 
scalar held is given by the relation, 


where p is the redefined canonical 
in the Einstein frame, 


dp 

d(j) 


\ 


12 2 + 



12 2 


scalar held. With these techniques we 


( 2 . 11 ) 

finally obtain the Lagrangian 
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where the superscript ( E ) represents the Einstein frame. The effective potential (12.131) has the £- and 
/^-dependences in addition to the A-dependence. The first one comes from the Weyl transformation and 
the second one from the quantum corrections. 


3 Inflationary parameters 


CMB fluctuations arise from the quantum fluctuations of the inhaton held. The measurements of CMB 
huctuations directly restrict the inhationary parameters. Within the context of the slow-roll scenario, the 
inhationary parameters are fully represented by means of the inhaton potential. The e-folding number, 
N, the slow-roll parameters, e and 77 , are formulated as [7], 
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where the upper and lower limits of the integral, (j>N and 4 > en d, are the field configurations when the 
slow-rolls scenario starts and brakes respectively. It is assumed that the inflaton slowly rolls downhill 
of the potential from V(^>n) to V{(j) en d)- The formula of scalar (or density) fluctuations <5 is given in 
Ref. US]- We modify it including the field redefinition (12.111) . 
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where the constant, C, appearing in Ref. [23] is taken to be 1. 

We apply these formulae to the two-loop effective potential. Substituting the effective potential (12.131) 
into Eqs. (EED-flUD, the inflationary parameters in the two-loop q l 4 model read 
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Under the slow-roll approximation the spectral index, n s , and the tensor-to-scalar ratio, r, are given by 


n s — 1 + 2 6 e|^_^ , (3-12) 

r = 16 e| 0=0JV - (3.13) 

These expressions are valid only if the slow-roll parameters are small enough, £< 1 and |?y| <C 1. The 
field configuration <j> en d is fixed by the boundary condition where the slow-roll parameters exceed the 
value of order one. 

It is instructive to take the small coupling limit, Ao —> 0. At this limit we can analytically carry out 
the integration in Eq (13.51) and obtain the simpler forms for the inflationary parameters 
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The equations <13.14B - H3.1GI> coincide with the expressions in Ref. [6]. Apart from the scalar fluctuations 
(13.171) . the inflationary parameters are found to be independent of Ao at the limit, Ao —> 0. Thus 
the spectral index, n s , and the tensor-to-scalar ratio, r, are determined by setting the scalar-curvature 
coupling, £o- For 4>n 4>end, the slow-roll parameters reduce to 
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at the large £o limit. Thus the inflationary parameters, r and n s , approach to 
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This is just the universal attractor mentioned in Ref. [7]. On the other hand, the slow-roll parameters 
are also simplified at the small £o limit, 
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At this limit the asymptotic values of r and n s are given by 
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As is pointed out in Ref. [3D], Eq. (13.211) does not represent the attractor. It is the model-dependent 
result in the various chaotic inflation models with the minimal curvature coupling. 


4 Lyapunov stability 


The spacetime of the inflationary universe is thought to be approximately de Sitter space. The exponential 
evolution of the spacetime needs to become unstable in order to get exit to radiation/matter dominant 
eras. Thus to check the stability of the de Sitter solution is important. If we suppose the spatially flat 
Friedmann-Lemaitre-Robertson-Walker (FLRW) universe during the inflationary era, Lyapunov stability 
of the de Sitter solution can be evaluated by the following formula mmm, 
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where the function F(</>) is the inflaton potential and U((f>) is the term multiplied by the curvature R in 
the Lagrangian in Jordan frame. Notice that, in this paper, the Lyapunov stability Kf is evaluated in 
Jordan frame. The function V' and V" denote the first and second derivatives of V with respect to <fi. 
The de Sitter solution is unstable for a positive Kf and stable for a negative Kf. The fixed point, </>/, is 
determined by 
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Substituting the Weyl factor (12.101) to U and the effective potential (12.71) to V , we calculate Kf for the 
model (EU). There is no solution of Eq. (1421) for a positive scalar-curvature coupling £o- For a negative 
£o some held configurations satisfy Eq. m • The phase structure of Lyapunov stability for a negative £o 
is illustrated in Fig. [T] It is found that the de Sitter solution is unstable when the absolute value of the 
scalar-curvature coupling, £o, or the quartic scalar coupling, Ao, is small enough. The figure Q] is plotted 
for fi 2, = 1. The renormalization scale dependence is very small. The difference in the results is only a 
few percent for a wide range of scales from fi 2 = 10 -6 to 10 3 . 

It should be noted that the exit from an approximately de Sitter space can be achieved through the 
reheating process. At this era we should carefully consider the interactions of the inflaton to light particles 
which assure the reheating of the universe [32 . The field configuration, <j>f, should have a value between 
<f>N and (j) en d as long as one assumes the slow-roll inflation scenario. In the considered model we cannot 
get the appropriate values of 4>n and <t> en d for a negative £o- It is unidentified whether the solution, <f>f, 
lies in the inflationary era or not. However, we believe that the above result shows something true and 
the de Sitter solution is unstable for a positive £o- In the next section we consider only positive values of 
the scalar-curvature coupling, £o> and numerically evaluate the inflationary parameters. 


5 

















-100 
-200 
-300 
^5 -400 
-500 
-600 
-700 

0.0 0.2 0.4 0.6 0.8 1.0 

Ao 


- Unstable 
: (*/><>) 

Stable 
[K f < 0) 



Figure 1: The phase structure for the Lyapunov stability for /i 2 = 1. 


5 Numerical results 


In Sec.3 we have analytically derived the explicit expressions for the inflationary parameters. These 
parameters depend on the scalar-curvature coupling, £o, the scalar quartic coupling, Ao, and the renor¬ 
malization scale, /i. These parameters also get the dependence of the e-folding number N through cfij^. 
In this section we numerically investigate the behavior of the spectral index, n s , and the tensor-to-scalar 
ratio, r, with varying the four model parameters, £o, Ao, /i 2 and N. The results are compared with the 
ones at the leading-log level obtained in Ref. [2]. 

First we determine the field configurations at the beginning and the end of the slow-roll scenario, (fix 
and (fiend, • At the end of the scenario the orders of the slow-roll parameters exceed unity. Thus the end 
of the scenario can be determined by 


I (j) — <j>end 


(5.1) 


Numerically performing the integration in Eq. (ESI) with N fixed, we obtain cfiisr. The inflationary 
parameters should be evaluated at the horizon crossing, <fi = cfipf, and so get the N dependence. In Fig. [2] 
we show the behavior of the field configuration, <fi^, as a function of the scalar-curvature coupling, £o- 
The field configuration, (fijy, strongly depends on £o> and is much greater than <fi en d, as seen in Fig. [2] 
The solid and dashed lines represent the results for the two-loop potential and leading-log potential, 
respectively. The difference between the results from the leading-log potential and the two-loop potential 
is very small, less than 2%. 

The slow-roll parameters, e\^ =(t)N and r )are numerically calculated by Eqs. (13.61) and (13.71) . 
Then we obtain the spectral index, n a , and tensor-to-scalar ratio, r from Eqs. (13.121) and (13.131) . In 
Fig. [3] we illustrate the ^-dependence of n s and r for 0.005 < £o < 10 5 . It is observed that the two-loop 
correction to n s is larger than the leading-log correction. The difference increases for a larger Ao- It is 
found that, whether for the leading-log or for the two-loop, the tensor-to-scalar ratio approaches r ~ 0.003 
independent of Ao, as £o increases. We cannot reproduce the attractor behaviors in HE] in the range, 
0.005 < £o < 10 5 . However it is not adequate to deny the existence of the attractors out of this range, 
i.e. a larger £ 0 10 5 . 

The renormalization scale dependence is illustrated in FigE We note that the lower value, /i 2 = 10“ 6 , 
corresponds to a typical GUT scale. It is observed that the tensor-to-scalar ratio, r, decreases as /i 2 
increases. This feature for the renormalization scale dependence is similar to that for £o i n both two-loop 
and leading-log results. For £o = 0.1 the tensor-to-scalar ratio takes the minimum value, r ~ 0.00874, 
at /x 2 = 10 3 in both results. It is worth noting that for a weak scalar-curvature coupling, £o = 0.1, 
n s increases as fj, 2 decreases in the two-loop result, while it is almost fixed for 10~ 6 < fj, 2 < 10 3 in the 
leading-log result. An opposite /independence of n s is observed for £ 0 = 10 5 in the two-loop result, unlike 
the leading-log result. It is considered that there is an optimized scale which minimizes the quantum 
corrections at a certain value of the inflaton field [TU]. It is normally set at the horizon exit, fi 2 = <fi 2 N , 
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Figure 2: Behavior of (j>N as a function of £o for N = 60 and fi 2 = 1. The solid and dashed lines are 
drawn for the two-loop and leading-log results, respectively. 
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Figure 3: Behavior of n s and r as a function of the scalar-curvature coupling £o hr the interval from 0.005 
to 10 5 for N = 60, /i 2 = 1 and Ao = 0.05,1. The solid and dashed lines are drawn for the two-loop and 
leading-log results respectively. 
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Figure 4: Behavior of n s and r as a function of the renormalization scale, /z 2 , in the interval from 10~ 6 
to 10 3 for N = 60, Ao = 0.5 and Co = 0.1,10 5 . The solid and dashed lines are drawn for the two-loop and 
leading-log results respectively. 


or at the end of the slow-roll scenario, /z 2 = (j) 2 N [53]. These optimized scales are also plotted in FigQ] A 
smaller two-loop contribution is observed for the renormalization scale, /z 2 = (jr N . 

Fig. [5] represents the Ao-dependence of the inflationary parameters. Since the two-loop effective 
potential (12.311 is perturbatively calculated, we vary Ao in the interval from 0 to 1. It is found that the 
tensor-to-scalar ratio, r, increases with the growth of Ao- The Ao-dependence of n s is enhanced for a 
larger Co in the two-loop result, which is again in contrast to the leading-log result. 

We plot the e-folding number dependence in Fig. [6] shifting N from 50 to 60, in which it is sufficient 
to solve the flatness or horizon problem. We see the tensor-to-scalar ratio, r, decreases as the e-folding 
number, N, increases. The range of r is almost equivalent for both two-loop and leading-log results. 
Furthermore we verify the suppression of 7V-dependence for a larger Co- 

The inflationary parameters are constrained from the observations. Current observational constraints 
on the inflationary parameters are given in Table [1] The scalar fluctuations and the spectral index have 
been observed by Planck and a recent joint analysis of BICEP2/Keck with Planck. Only an upper limit 
is obtained for the tensor-to-scalar ratio, r. 


Observations 

5 x 10 5 

n s 

r 

Planck + WP + lensing 

4.93 

0.9653 ±0.0069 

< 0.13 

BICEP2/Keck + Planck 

— 

— 

< 0.12 


Table 1: Constraints on inflationary parameters by Planck BJ 1311 and BICEP2/Keck with Planck [5]. 


We try to explain these constraints under the two loop effective potential. For the intervals of the 
four model parameters 

0.005 < Co < 10 5 , 0 < A 0 < 1, 10" 6 < n 2 < 10 3 , 50 < N < 60, (5.2) 

we search suitable parameter sets. As for the effective scalar-dependent quartic coupling we also keep 

0 < A < 1. (5.3) 

It is possible to reproduce the values of the scalar fluctuations and spectral index obtained by Planck 
and WMAP. In Table [2] we summarize the parameter sets which are consistent with the observational 
data for the scalar fluctuations and the spectral index. The tensor-to-scalar ratio does not exceed the 
upper limit by BICEP2/Keck and Planck. It is found that the scalar-curvature coupling, Co, and the 
renormalization scale, /z 2 , play a crucial role in determining the inflationary parameters. It is noted again 
that the case with /z 2 = 10~ 6 is appropriate for a typical energy scale of inflation, about 10 15 GeV. 
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Figure 5: Behavior of n s and r as a function of the 
scalar quartic coupling, Ao, in the interval from 
0.05 to 1 for TV = 60, fi 2 = 1 and £ 0 = 0.1,10 5 . 
The solid and dashed lines are drawn for the two- 
loop and leading-log results respectively. 


Figure 6: Behavior of n s and rasa function of the 
e-folding number, TV, in the interval from 50 to 60 
for /i 2 = 1, Ao = 0.5 and £o = 0.1,10 5 . The solid 
and dashed lines are drawn for the two-loop and 
leading-log results respectively. 


(Co, A 0 ,^ 2 ,TV) 

A 

5 x 10 5 

n s 

r 

(44.4,0.0383,10- e , 60) 
(44.4,0.0385,10“ 6 ,50) 

0.03850 

0.03870 

4.9332 

4.9282 

0.96529 

0.96531 

0.002991 

0.004227 

(91.8,0.17,1,60) 

(93.6,0.1775,1,50) 

0.16998 

0.17742 

4.9306 

4.9329 

0.96534 

0.96531 

0.003038 

0.004285 

(150.6,0.5,13.9,60) 

(151.5,0.5,11.65,50) 

0.49264 

0.49262 

4.9307 

4.9300 

0.96529 

0.96529 

0.003166 

0.004433 

(200.6,1,27.38,60) 

(203.6,1,22.88,50) 

0.96309 

0.96293 

4.9306 

4.9302 

0.96533 

0.96534 

0.003350 

0.004649 


Table 2: Our results of inflationary parameters. 
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6 Conclusions 


In this paper we have studied a two-loop </> 4 theory with non-minimal scalar-curvature coupling. Assuming 
the slow-roll scenario, we have analytically derived the expressions for the inflationary parameters, the 
scalar fluctuations, 5, the spectral index, n s , and the tensor-to-scalar ratio, r. The Lyaounov stability 
has been evaluated for a negative £o- In this region the de Sitter solutions become unstable when the 
absolute values of the couplings decrease. 

The behavior for n s and r has been systematically evaluated as a function of the couplings, the 
renormalization scale and the e-folding number. It is found that the two-loop corrections on n s are larger 
than the leading-log corrections. Besides, it is demonstrated that a larger Ao- and /^-dependence of n s 
is observed as the scalar-curvature coupling, £q, increases. We have also observed a similar feature for 
£o-, fi 2 - and IV-dependence of r. Concretely the tensor-to-scalar ratio r is suppressed with the growth of 
£o, I* 2 and N. 

The loop corrections introduce an unavoidable dependence on the renormalization scale. It is observed 
that the two-loop contribution is softened for p 2 = (j) 2 N . It should be noted that we can tune the couplings, 
Ao and £o consistent with the current observational constraints for n s and r. In SM Higgs inflation the 
quartic coupling vanishes near the Planck scale m- a small positive value of the non-minimal scalar 
curvature coupling is favored for a simple scalar quadratic potential [35] . 

It has been reported in Ref. [7] that there is an attractor in the (n s ,r) plane at the large £o limit. 
We find that such attractor behavior does not appear in our model until £o = 10 5 6 . But the results seem 
to approach the attractor. There is a possibility that a larger £ 0 is necessary to obviously observe the 
convergence to the attractor due to the quantum effects. It has been pointed that multi-field models with 
quartic self-couplings quickly relax to the single-field attractor for a wide range of couplings and initial 
conditions [36]. Thus it is interesting to extend our approach to multi-field models. 

It is mentioned that the formulation of scalar-tensor theories and the theoretical predictions depend 
on the choice of the frames, Jordan or Einstein frame m ■ Note that calculation of inflationary param¬ 
eters and Lyapunov stability may look to depend on the Einstein or Jordan frame under consideration. 
However, it has been proved in I38E2I that the obtained inflationary parameters are just the same in 
both frames. In Ref. [40] the induced inflation models are revealed to have another attractor at the 
small £ 0 limit. The radiative corrections may affect such a double attractor behavior. Note finally that 
using same approach one can study non-minimal two-loop inflation for Standard Model. This will be 
considered elsewhere. 
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